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ABSTRACT 

We investigate a magnetized plasma in which injected high energy gamma- 
rays annihilate on a soft photon field, that is provided by the synchrotron radi- 
ation of the created pairs. For a very wide range of magnetic fields, this process 
involves gamma-rays between 0.3 GeV and 30 TeV. We derive a simple dynam- 
ical system for this process, analyze its stability to runaway production of soft 
photons and paris, and find conditions for it to automatically quench by reaching 
a steady state with an optical depth to photon-photon annihilation larger than 
unity. We discuss applications to broad-band 7-ray emitters, in particular su- 
permassive black holes. Automatic quenching limits the gamma-ray luminosity 
of these objects and predicts substantial pair loading of the jets of less active 



sources. 



Subject headings: radiation mechanisms: non-thermal 
galaxies: active 



1. Introduction 



gamma-rays: theory 



Astrophysical sources of high energy X- ray and 7-ray radiation are often comp a ct enough 
for significant photon-photon annihilation (jJelleyl Il966l : iBonometto fc Reed Il97ll : iHerterich 



1974 



Guilbert et al.l Il983l ; iKazanad Il984l ). leading to electromagnetic cascades. In many 



situations, the cascades can be assumed to be linear, in the sense that the density of target 
photons is fixe d. This is particularly convenient for a Monte-Carlo computation of the 
spectrum (e.g., lAharonian et al.lll985t iProtherod Il986l ) . Nonlinearity arises when photons 
produced by the cascading electron-positron pairs are themselves targets in the photon- 
photon pair-production process. This has been considered mostly in connection with compact 
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X-ray sources (jSvenssonl 1 198 71 ) . where photons that are inverse- Compton scattered by the 
nonthermal pairs form the target population. 



Using simplified kinetic equations, IZdziarski fc Lightmanl (119851 ) studied the steady state 
radiative spectra of non-thermal electrons in active galactic nuclei subjected to inverse- 
Compton and synchrotron energy losses, taking photon-photon annihilation and the result- 
ing production of a first generation of secondary electron-positron p airs into account self- 
consis tently. A similar discussion was presented by lFabian et al.l ( 19861 ) and lLightman fc Zdziarski 



( 119871 ). who investigated the spectra of such sources for a wide range of parameters and 
included the effects of thermal pairs. Facilitated by the dramatic improvement in comput- 
ing resources, algorithms for the computation of time-dependent solutions of the kinetic 
equations describing these proc e sses have been d eveloped by several groups (e.g., ICoppi 



19921 : iMastichiadis fc Kirkl Il995 [_| stern et all 1 19951) an d are routinely used in source mod 
eling ( Konopelko et al. Il999l ; Krawczynski et al. 2002 ) . Recently, however , atten tion has 
focused on intrinsically nonlinear effects. For example, iKirk fc Mastichiadisl (119921 ) studied 
the conditions under which ultrarelativistic protons initiate runaway pair-production. For 
active galactic nuclei, t he parameters are ex t reme, but gamma-ray bursts pr ovide a more 
pro mising environment (IKazanas et al.l 120021 ; IMastichiadis fc Kazanad 120061 ) . ISternl (120031 ) 



and Stern Sz Poutanen (l2006h have found a strong feedback mechanism that may mediate 



relativistic shocks, or, in relativistic jets, breed high energy 7-rays. 

In this Letter we analyze the conditions under which a magnetized system, with an 
embedded source of photons of energy e = em e c 2 ^> 1 can sustain pair production au- 
tonomously. The nonlinear process that makes this possible is annihilation on soft photons 
with energies eo = eo m e c 2 <C 1, produced as the synchrotron emission of the created electron- 
positron pairs. To describe the dynamics of this process, we simplify the description to a 
system of three, first-order, ordinary differential equations for the evolution on timescales 
larger than the light-crossing time of the source, of the populations of hard photons, soft 
photons and electrons. To this end, we assume an isotropic distribution of the photon fields 
and relativistic pairs in the source rest frame. The described process operates when the 
magnetic field B is such that 



B 



6 B„ e 



where B CT = m 2 e c 3 /he « 4.4x 10 13 G. This implies B = 265 (e/10 6 ) jiG = 35 (e/TeV)~ d //G. 
Such systems are important in high energy astrophysics, because for an extremely wide 
range of magnetic field strengths (B = 10~ 9 G— 10 6 G) Eq. ([T]) is satisfied for photon ener- 
gies between e ~ 0.3 GeV and 30TeV, that are accessible to space or ground based 7-ray 
observatories. 



The reduction of a similar complex system of relativistic particles and photons to a 
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dy namical system of two o rdinary differential equations was discussed in an early paper 
by IStern &: Svenssonl (Il99ll ). who found limit-cycle solutions. In principle, the dynamical 
system we derive is of sufficiently high order (third) to display much richer "chaotic" behavior. 
However, the physical processes we discuss are strongly dissipative; in the system studied, 
we find no physically acceptable oscillatory or quasi-oscillatory solutions and no chaotic 
regime. Nevertheless, the system has interesting stable solutions with implications for the 
7-ray spectrum of supermassive black holes, and the electron-positron pair loading of their 
jets. 



2. Equations 



We restrict our analysis to the high energy and soft photon populations that interact 
with ultrarelativistic pairs (henceforth "electrons") of characteristic Lorentz factor 7. In the 
systems of interest, the radiative cooling time scale of these electrons is much shorter than 
the light-crossing time of the emitting volume, and we assume that no other process such as 
acceleration, or interaction with background particles or waves competes with the cooling. 
Although our formulation does not enable us to extract detailed spectral information, we 
discuss the implications for the case of a broad-band spectrum of high energy photons. 

We denote the total number of high energy photons, soft photons, and the created 
electrons in our source volume V by N, N Q and N e , and define the effective logarithmic 
number densities asJV(e) = N/eV, Ao(eo) = N /eoV, A4(7) = N C /^V. The kinetic equations 
are then 



dAT(e) 

dt 
dA/- (e ) 

dt 
dA4(7) 
dt 



^o esc (eo)+A/ syn ( eo ) , 
AC 7 (7)+A/ syn (7)+A/ ic (7) 



(2) 



sec 



Zdziarski fc Lightmanl Il985l ; iFabian et al.l 119861 ) . The terms on the right-hand sides 
of the Eq. (j2J) are: (i) the injection rate Q(e) = L in j/(e 2 m c c 2 V) of high energy pho- 
tons into the emission volume V = 7ri? 3 , where is the high energy luminosity intrin- 
sic to the source. The annihilation rate of these photons A/" 77 (e) = — ca 77 A/"(e), with 
a 77 = (2/3) o"t e _1 A/o(e ), as follows from the approximate exp ression for the photon-photo n 



annihilation cross section, a. 
and their escape rate J\ ' 



77 

esc/ , 



= (l/3)(T T e 5[e - (2/e)] flZdziarski fc Lightmanlll985f ). 
-J\f(e)/t esc with t esc = R/c (this defines the effective 



radius R). (ii) The soft photon production rate by synchrotron emission eoA/Q iyn (eo) 
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7-A/" e (7) |7l 

„2 



syn 



m e c 2 /(87r). the syn- 



47r j t0 /(m e c 2 ), with the synchrotron emissivity eoj eo 

chrotron energy losses term |7| syn = (4 c<Tt /3 m e c 2 ) Ub 7 2 , and the magnetic field energy 
density [7b = B 2 /Stc, and their escape rate A/"n sc (en) = — A/"n(en)/t PS f (in) The pair injection 
rate A/^ 7 (7) = -47v^ 7 (e) (since 7 » e/2, see ICoppi fc Blandfordl Il99ph . Their radiative 
cooling rate, by synchrotron emission, A/" e syn (7) = — Af e (j)/tsyn where t syn = 7/|7| sy n, and 
by inverse-Compton scattering A/" e lc (7) = — Ae(7)/^i c with i ic = T / 1 "T I ic - both modeled here 
as catastrophic escape from the analyzed energy range. The inverse-Compton loss term is 
|7| ic = (4 co"t/3 m e c 2 ) C/ 7 2 ?7~ 1 , where U = m e c 2 Af (e ) is the energy density of the soft 

1 and 
The influence 



photons and the factor rj « 5 1 ' 5 takes account of Kl ein-Nishina effects, sin ce 760 
we do not consider the detailed spectral shape 



sec 



Moderski et all 120051 ) 



of inverse-Compton scattering and annihilation on the soft photon field is negligibly small 
provided 7 ^> 1 and e 1. 



or B), R and L in j, and three variables, 



Our system thus has three free parameters, e 
AA(e), A/o(eo) and M e (l)- It is convenient to express the parameters in terms of dimensionless 
quantities: the high-energy injection and the magnetic compactnesses, 



e^inj 0t 



-inj 



ATI m e c 3 R 



and 



m P cr 



(3) 



Note that these qua ntities contain a factor e not present in the conventional definitions 
(IGuilbert et al.lll983l ). Instead of the soft photon number density we use the optical depth 
for photon-photon annihilation n = t 77 = 2 cr T RJ\f (e ) / (3e), and instead of the hard photon 
number density, we use the ratio of the hard photon energy density to the energy density 



of soft photons needed to give r 77 = 1: n = UJ Uq\ t 



e 3 <TT-RA/'(e)/6. Similarly, for 



the electron density we use the ratio of the energy density to that of soft photons required 
for r 77 = 1: n c = U c / Uq\ t =1 = e 3 <Tt -RA/" e (7)/24, and measure time in units of the light 
crossing time r = ct/R. This leads to the autonomous dynamical system: 



dn 

d7 

dn 

dr 

dn e 

dr 



-nn - n + - £ inj 



2 
3 



-n + - 4s n e 

n n - - £b n e n e n 

6 7] 



(4) 



3. Stationary solutions and their stability 



For all £inj > and Ib > there exists a trivial stationary solution of the system (j4]): 
n = 2£ in j/3, n = 0, n e = 0, corresponding to free propagation of the hard photons through 
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a medium in which pairs and soft photons are completely absent. The eigenvalues of the 
Jacobian matrix evaluated for this solution are all real. For £ in j < 3 they are all negative, so 
that the solution is stable in this region. This means that the nonlinearity is not sufficient 
for the system to initiate and sustain a population of pairs if £- in ^ < 3. Furthermore, in 
this region, no other physically acceptable (no > 0, n c > 0) stationary solution exists. For 
£i n j > 3, at least one of the eigenvalues has a positive real part, so that the solution in which 
pairs are absent is unstable. However, in this region a second stationary solution emerges, 

with n = 2£ in j/ [3 (1 + no)], n c = 3n /£ B and n = 



-6 - rje B + J (r]£ B - 6) 2 + 8 V £ B £ ini 



/12. 



Using the Routh-Hurwitz conditions (IGradshteyn fc Ryzhik! Il980l . 17.715), one can show 
that the real parts of the eigenvalues of the Jacobian matrix corresponding to this physical 
solution are negative, so that the solution is stable for all values of £ B > 0. All physically 
acceptable solutions appear to approach this stable one. We find no chaotic regime. We also 
have no physical grounds on which to expect such behavior within the current model. 

Figure [JJ illustrates the properties of the stable solutions of this system on the £ B -£m) 
plane. The horizontal dotted line £- m j = 3 bounds from above the region (shaded dark gray) 
in which the stable solution contains neither pairs nor soft photons, so that r 77 = 0. The 
solid line, on which £ in j = 6 + 36/(r/£ B ), is the locus of points at which the stable solution has 
unit optical depth to absorption of high energy photons: n Q = r 77 = 1. Strictly speaking, 
our treatment of the escape probability of hard photons is valid only in the optically thin 
case r 77 < 1, which is shaded light gray. However, for £ in j > 6 + 36/(r]£ B ), (white area) 
we expect on physical grounds that the system tends spontaneously to an optically thick 
state with almost total absorption of the hard photons. Stable solutions which have equal 
energy densities of soft and hard photons (n = no) occur formally when £ in j = 3r]£ B (3r]£ B — 
12)/(r]£ B — 12) 2 , shown as a thick dashed line. This line lies entirely within the optically 
thick region. Thus, our system is valid only for solutions in which the energy density of soft 
photons is small compared to that in hard photons. Finally, stable solutions in which the 
electron and magnetic field energy densities are equal occur when U C /U B = 6n e /£ B = 1, 
shown as a thick dotted line. Solutions which are particle dominated lie above this line in 
the hatched region of Fig. [TJ It is interesting to note that when the optical depth for photon- 
photon annihilation is less than unity, r 77 < 1, stable solutions can be found that are pair 
dominated, in the sense that the electron energy density exceeds that in the magnetic field. 
These solutions occur in the intersection of the hatched and light gray areas of Fig. [TJ They 
imply efficient pair- loading of the system and are particularly relevant for Ib C 1, where 
they occur for all relevant injection compactnesses £ in ] > 3. 
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4. Applications 

The electrons created in the photon-photon annihilation process cool due to synchrotron 
and inverse- Co mpton energy losses. In the analysis presented so far, where a monoenergetic 
high energy photon field was considered, with fixed photon energy e = e*, these losses were 
modeled as the escape of particles from the relevant energy range. In a more precise treat- 
ment, the electrons do not escape catastrophically, but cool (due to synchrotron emission and 
inverse Compton scattering). Consequently, below the injection energy 7* = e*/2 the steady 
state spectrum of the cooled electrons i s N^jj < 7*) oc 7~ 2 , leading to a synchrotron spec- 



trum at photon energies e < = 2/e* (jKardashevlll962l ) . This corresponds to an emissivity 
j eo <e^ oc eo -1 / 2 . In turn these photons absorb high-energy radiation with photon energies 
e > e*. Since the optical depth for the photon-photon annihilation is r 77 oc e Ao(eo) oc j eo , 
the condition for r 77 > 1 at e = e* ensures automatically r 77 > 1 for any e > e*. In 
other words, if the parameters of the considered system are such that it becomes optically 
thick for the high energy photons with e*, then all the emission at higher photon energies 
e > e* will also be absorbed. For an intrinsic luminosity L in j corresponding to the photon 
energies e*, the appropriate condition for r 77 (e > e*) > 1 can be therefore expressed as 
Anj > £ cri t = (lO- 6 ^ eV 4 + 4.4 x 10 5 i? pc /4ev) x 10 40 erg/s . 

Supermassive black holes (SMBHs) in active galactic nuclei are an interesting appli- 
catio n of the process d i scussed here. These objects may accelerate protons to high ener- 



gies (ISikora etal.1 119871 ; iBoldt &: Ghoshlll999l ). generating very high energy 7-ray emission 
(jLevinsonll2000l ) . The magnetic field intensity in such systems is high: assuming equipartition 
between magnetic and accretion-produced radiation gives B = (i?/r g ) _1 (Lacc/^Edd) 1 ^ 2 (-^Edd/^gc) 
where r g = GAi/c 2 pa 1.5 x 10 13 A^scm is the Schwarzchild radius of the black hole 
with mass M. = .Ms 10 8 M , L acc is the accretion luminosity produced within radius R, 
and L Edd = An G M.m v c/ k 1.3 x 10 46 Ai 8 ergs -1 is the corresponding Eddington lu- 
minosity. SMBHs in AGNs typically have masses Ai$ ~ (10~ 2 — 10 2 ) and luminosities 
A = L acc /LEdd ~ (10 -7 — 1). However, lower-mass black holes in the centers of 'normal' 
galaxies, with Ai$ < 10~ 2 , also exhibit nuclear activity with very low accretion luminosi- 
ties A ~ 10~ 10 (e.g., Sgr A* in the Galactic center). Assuming that the magnetic field 
scaling holds for all these objects within R ~ few x r g , one finds B ~ 10 4 (A/A^ 8 ) 1/2 G, 
and, hence [Lc^/lO^ergs- 1 ] ~ 3 x 1(T 6 (A/A-l 8 )^ 2/3 + 5 M 8 {X/M 8 ) 1/6 at photon energies 



1/2 



(\/M 8 y 1/6 GeV. 



The critical luminosity and the corresponding values of e* are plotted in the upper (left- 
hand ordinate) and lower panels of Figure [2J respectively, for different black hole masses. 
In the case of galactic nuclei hosting the most massive black holes Ai 8 > 1 no steady 
(on the light-crossing timescale) high energy 7-ray emission at photon energies e > e* ~ 
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(1 - 100) GeV can emerge if [e*L £ *] > (10 41 - 10 42 )erg s _1 . For M 8 < 1, L C rit is much 
lower, ~ (10 38 — 10 39 ) erg s _1 . This does not mean that lower luminosity sources are able to 
emit high-energy gamma-rays, since these may also be absorbed on the 'external' low-energy 
radiation of the accreting matter. The relevant optical depth is r ext = (1/3) ctt RU ex t/e ext , 
where U ex t is the monochromatic energy density of the accretion-related emission at e cx t = 
2mgC 4 /e*. A conservative estimate is U ext = ( L ajCC /(4:7i R 2 c), with ( »s 0.1, since, for 
e* ~ (1 — 100) GeV, it is the low-energy photons e ext ~ (0.01 — l)keV that are relevant. 
Using Eq. ([T]) and the assumed scaling of the magnetic field r cxt « 2.7 x 10 4 A (A/A^s) - ^ 6 - 
The upper panel of Fig. [2] shows r ex t (thin lines, right-hand ordinate), and one sees that 
automatic quenching of high-energy 7-ray emission generated in the immediate vicinities of 
SMBHs competes with the opacity due to emission of the accreting matter only in the case 
of low accretion rates, A < 10 -5 . 

The loading of the system with electron-positron pairs is another interesting aspect of 
automatic quenching. Defining Aj n j = L in j/Z/Edd, and again assuming magnetic and radiation 
pressures are equal, the compactnesses (Eq. [3]) become £^ ~ 10 6 Aj n j (\/ M.%)~ 1 ^ and £b ~ 
10 6 A (\/M 8 )~ 1/s . Fig. Q] shows that the energy density of the pairs substantially exceeds 
that in the magnetic field if £ in j > 3, £b < 1, and r ext < 1. In the range of A and A4 8 
of interest for SMBHs one finds (X/M 8 ) /D ~ (1 - 100). Consequently, pair loading is 
important for those lower luminosity sources, whose gamma-ray output is comparable to, or 
less than, the accretion luminosity A < (10~ 8 — 10~ 6 ) < A in j. 



5. Conclusions 

The resonance-like behavior of the photon-photon annihilation cross section means that, 
for magnetic fields between B = 10~ 9 G and 10 6 G, 7-ray photons between e ~ 0.3 GeV and 
30 TeV (see Eq. [T|) can be absorbed on the synchrotron radiation of the created pairs. This 
energy range is now, or will shortly, be accessible to the AGILE and GLAST satellites, 
together with the Cherenkov Telescopes HESS, MAGIC, CANGAROO, and VERITAS The 
nonlinear dynamics of these 7-ray sources render them unstable to runaway pair production 
when the injection compactness £i n j > 3 (see Eq. |3]). The spontaneously produced optical 
depth to absorption saturates at value greater than unity for photons with energies above 
£, if Anj > 6 + 36/(77 £b); where 77 5 1 ' 5 and £b is the magnetic compactness (see Eq. [3]). 
In this case the source automatically quenches its own 7-ray emission. As a corollary, we 
find that systems with small £b < 1 but large £- m j > 3 suffer significant pair loading. This 
implies that 7-ray sources in the vicinity of low accretion-rate supermassive black holes can 
be automatically quenched and their jets loaded with electron-positron pairs. 
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Fig. 1. — Properties of the stable solutions to the dynamical system of Eq. as functions of 
the injection and magnetic compactnesses, Eq. The dark gray shading denotes solutions 
with no electrons or soft photons, light gray shading denotes solutions with optical depth to 
photon-photon absorption less than unity and in the white area strong absorption occurs. 
The hatched area denotes solutions in which the energy density of the created pairs exceeds 
that in the magnetic field. The dashed line is the locus of points at which the energy density 
in soft photons formally reaches equipartition with that in hard photons, the dotted line 
that of equipartition between pairs and magnetic field. 
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Fig. 2. — Upper panel: the critical injection luminosity L crit that automatically generates a 
pair density and soft photon field of unit optical depth to photon-photon absorption (thick 
lines, left-hand ordinate), and the opacity of e* photons emitted by accreting matter (thin 
lines, right-hand ordinate), both as functions of the accretion luminosity L acc /L Edd . Lower 
panel: the corresponding maximum gamma-ray energy e*. In each case four lines are plotted, 
corresponding to different masses of the SMBH. 



